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Introduction 

The discovery of paraffin based fuels has lead to renewed interest in hybrid rocket 
research. Experiments have shown that they bum 3-5 times faster than conventional 
hybrid fuels. High thrust level that would have required a multi-port design in the past 
can now be achieved with a single-port motor. While tests performed in Stanford and 
NASA Ames have demonstrated the paraffin hybrids to be a promising technology, one 
of the major challenges has been the relatively low efficiency. The c* efficiency has 
ranged between 80% and 90% in experiments conducted at the Ames Hybrid Combustion 
Facility (HCF). The test motor in these experiments had a 45 inch long fuel grain with the 
initial port diameter ranging between 3 and 5 _ inches. The c* efficiency is defined as the 
ratio of measured and theoretical characteristic velocities and is related to how 
completely the fuel and oxidizer are converted to combustion products. A low efficiency 
means that the reactants bum incompletely, and the reaction does not release the 
maximum possible amount of energy. 

Two species react incompletely if either their chemical reaction rate is slow or if they are 
not mixed adequately. Chemical reactions take place through molecular collisions, and if 
the two species do not come into contact because of insufficient mixing, they can not 
form products. Typical fuels and oxidizers exhibit fast reaction rates. The characteristic 
time scale of chemical reactions is much smaller than the time scales related to 
convection and diffusion. In other words, the Damk hler number, Da=_ P h ys /_chem, where 
_phys and _ C hem are characteristic physical and chemical time scales, is large. In this limit 
of large Da, chemical kinetics is not thought to influence the efficiency. Therefore, a high 
efficiency can be accomplished by effective mixing. It can be challenging to achieve 
good mixing between the reactants in a hybrid motor. Fuel and oxidizer are stored 
separately and in different phases. The solid fuel must vaporize and the reactants must 
diffuse to meet each other before chemical reactions can take place. In contrast, the fuel 
and oxidizer in a solid rocket are premixed, and in liquid-liquid systems, the reactants can 
be injected into the combustion chamber in a way that promotes mixing. Consequently, 
hybrids are generally less efficient than solid or liquid systems. It should be noted, 
however, that while storing the reactants separately makes it difficult achieve high 
efficiency, the same feature also makes hybrid rockets inherently safe which is one of 
their major advantages over other types of rockets. 

The efficiency of a hybrid rocket is generally found to increase with the length of the 
motor. The residence time increases with the length of the motor, allowing more time for 
the reactants to diffuse towards each other. However, simply making a motor longer to 
increase its efficiency is not a practical solution, and alternate approaches must be 


considered. The focus of the present research is to find a way to mix the fuel and oxidizer 
better that is practical and simple to implement in a hybrid motor. The approach taken is 
to alter the flow field by generating disturbances to the oxidizer flow upstream of the 
combustion chamber, and to look for a configuration that promotes mixing. Passive 
mixing devices inserted to the oxidizer flow path or changes in the pre-combustion 
chamber area geometry can be used to produce such disturbances. 

A cold flow experiment was designed for investigating the mixing process in a hybrid 
rocket. Operating a chemically inert experiment is simpler, safer and less costly than 
firing the real motor. While heat release will influence the flow field in the chemically 
reacting case, the flow patterns should be qualitatively similar to those in the cold flow 
situation, and the mixing aspect of the flow field can be captured in the non-reacting 
environment. A configuration that produces a high degree of mixing in the cold flow case 
is expected to do so in the presence of chemical reactions. 

This document describes the cold flow experiment performed to study mixing in a hybrid 
motor, and a numerical simulation of the experiment. 



Figure 1 : Mixing experiment apparatus 


Experiment setup 

Figure 1 shows a photograph of the experimental apparatus. Both the fuel and the 
oxidizer are simulated with air. The fluid representing oxidizer flows axially in a pipe, 
while the fluid representing fuel is injected from the pipe walls. The radial flow is 
maintained at a temperature slightly higher than that of the axial flow. Mixing is inferred 
from temperature profile measurements taken at the pipe exit. 

Figure 2 shows a schematic of the apparatus. The axial flow is driven by a high pressure 
blower, and the flow’ rate is controlled by adjusting the speed of a motor driving the 
blower. The fluid first flows through a long, 3 inch diameter inlet pipe. This will produce 
a fully developed velocity profile at the downstream end of the inlet pipe. The pipe also 
contains honeycomb and screens for flow conditioning. After the inlet pipe, the flow 
expands through what would be the pre-combustion chamber area in a hybrid motor, and 
finally enters the main pipe, which corresponds to the fuel port of a motor. The diameter 
of the pipe is 4 inches, and it is 44 inches long. The fluid representing the fuel comes 



from a compressed air supply. The flow rate is controlled by a choked orifice. Before 
being injected into the axial flow, the fluid enters an outer chamber surrounding the main 
pipe through a distribution tube. The tube contains holes that increase in size along its 
length to distribute the fluid into the chamber as uniformly as possible. The fluid then 
enters the main pipe through the wall. The pipe wall is a perforated metal cylinder, and it 
is surrounded by layers of foam and low-permeability cloth. The foam and the cloth 
create a significant pressure difference between the interior of the pipe and the outer 
chamber. The pressure drop between the pipe centerline and the outer chamber for a 
typical run is about 350 Pa, while the difference between the maximum and minimum 
pressure inside the pipe is of the order 20 Pa. This ensures that the radial flow at the wall 
is as uniform as possible along the length of the pipe. After entering the pipe, the axial 
and radial flows can interact and mix, and finally exit the pipe. 

The temperature of each air stream is actively controlled. A 1000 W heater located 
downstream of the choked orifice controls the temperature of the radially injected fluid. 
The heater is connected in series with a solid state relay. An Omega CN8501 temperature 
controller unit switches the relay on and off, maintaining the fluid at a desired setpoint. 
The temperature of the axial stream is set by another controller of the same model. A 
500W heater is located upstream of the blower. The axial flow is typically maintained at 
a degree or two above the ambient temperature, and the radial flow is heated by an 
additional 10 K. 

The geometry of the apparatus replicates that of the Ames Hybrid Combustion Facility, 
with the exception that it does not include a post-combustion chamber or a nozzle. The 
ratio of axial to radial mass flow rates in the experiment is maxi a i/m ra diai-2.4, 
corresponding to a typical oxidizer to fuel ratio in test firings at the Ames facility. The 
Reynolds number based on the pipe diameter and average velocity of the incoming axial 
flow is approximately 31 0 5 . 

The primary measurements have been mean temperature and velocity profiles at the exit 
of the pipe. Temperature profiles are recorded with an RTD probe, and velocity profiles 
are measured with a Pitot static probe used in conjunction with a high-precision Combist 
manometer. To measure the profiles, the appropriate sensor is mounted on a traverse. 
After each data point is recorded, a stepper motor moves the sensor by a distance set by 
an indexer. In addition to the movable sensor at the pipe exit, the experiment has four 
additional, permanently mounted RTD sensors. One of these sensors is mounted on the 
perforated wall of the pipe, and another one measures the temperature in the outer 
chamber. Two sensors monitor temperatures of the incoming air streams, and provide 
process values for the two temperature controllers. 



Figure 2: Schematic of experimental apparatus 



Computational Model 

The numerical simulation was a Reynolds Averaged Navier Stokes solution for the 
geometry of the experiment. Star-CD, a commercial CFD package was used for pre- 
processing, solving and post-processing. The turbulence model was the standard k-_ 
model with wall functions. The flow was assumed to be axisymmetric, and cylindrical 
coordinates were used in the model. The calculation was performed in a wedge extending 
10 degrees in the azimuthal direction. The solution was performed on two different grids. 
The base grid had approximately 12000 cells. This grid was refined by subdividing each 
cell into four cells, resulting in a grid with approximately 48000 cells. Views of a section 
of the base grid and refined grid are shown in figures 3 and 4. The computation required 
boundary conditions for both temperature and velocity. In the azimuthal direction, a 
symmetry boundary condition was specified for both variables on planes at +/- 5 degrees. 
A fixed axial velocity profile and a constant temperature were specified at the inlet, the 
temperature and radial velocity of the injected fluid were fixed at the pipe wall, and the 
mass flow rate was specified at the exit. 
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Figure 3: View of a part of the computational mesh, base grid 
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Figure 4: View of a part of the computational mesh, refined grid 


Results 

Comparison of Experimental and Computational Results 

Figure 5 shows a comparison of calculated and experimentally measured velocity and 
temperature profiles. The axial velocity at the pipe exit is plotted as a function of the 
radial distance from the pipe centerline. The temperature data is presented in the non- 
dimensional form, (T-T re f)/(T W aii-T re f). T re f is the temperature of the incoming axial flow, 
T wa n is the temperature of the perforated pipe wall, and T is the temperature measured at 
the pipe exit. The data corresponds to the baseline case with undisturbed axial flow. 
Computational results are shown for both the base grid and the refined grid. The match 
between the computational and experimental results is currently rather poor. Refining the 
grid did not bring the computational solution closer to the experimental results. This 
indicates that further refinement of the grid is not expected to produce a better match 
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between the two. Compared to the computation, the experimental results show a 
significantly larger difference between the centerline temperatures and the temperature at 
region near the wall. It appears that the computation is producing more turbulent 
diffusion than the experiment. Fluid injected from the wall is more likely to end up near 
the centerline according to the computation than according to the experimental data. 
Figure 6 shows a comparison between the calculated and experimentally measured static 
pressure as a function of axial distance in the pipe. This indicates that in the calculation, a 
larger pressure gradient is required to drive the flow through the pipe than what is 
observed in the experiment These discrepancies have not been satisfactorily explained. 
However, only the standard k-_ turbulence model has been attempted, and it is possible 
that another turbulence model might produce better results. It is also possible that the 
flow in question, particularly with the complex separation-reattachment zone at the 
expansion area upstream of the pipe is too complicated to be handled by a steady RANS 
calculation. More powerful computational tools may be required to improve the match 
between experiment and computation. 
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Figure 5: Comparison of experimental and computational results, velocity and 
temperature profiles at exit plane 
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Figure 6: Comparison of experimental and computational results, axial pressure profile 


Testing of Flow Disturbance Generators 

Figures 7-1 1 show images of the devices that have been tested. The devices are placed in 
the path of the axial stream, either at the downstream end of the inlet pipe or just 
upstream of the main pipe. The intent is to enhance mixing by generating disturbances to 
the axial flow. None of the devices tested so far have produced a significant deviation to 
the temperature profile compared to the baseline case, an example of which is 
demonstrated by the data in figure 1 2 showing the baseline temperature profile and the 
profile produced by one of the devices (shown in figure 7). 

One important point to be made is that a mean temperature measurement does not 
indicate whether or not the fluid is mixed at molecular level. The mean temperature of a 
flow of warm and cool fluid blobs may be the same as that of a flow of homogeneously 
mixed fluid. The mean temperature profile can only indicate how likely it is for the warm 
fluid entering through the wall to end up near the centerline, or for the cooler fluid to find 
its way from the centerline to a region near the wall. A measurement of the fine scale 
mixing should be incorporated to truly distinguish between mixed and unmixed fluid. 

The mean temperature profile by itself can only provide information about the bulk 
movement of the fluid in the radial direction. If the mean temperature between the wall 
and the centerline is small, the cool fluid has been transported towards the wall, and the 
warm fluid from the walls has moved towards the centerline. On the other hand, if the 
mean temperatures at the centerline and near the wall are very different, the cool fluid 



had stayed near the centerline, and the warm fluid still resides next to the wall. 
Ultimately, it might be desirable to achieve both an increase in the bulk transport of fluid 
in the radial direction as well as a high degree of mixing at the molecular level. However, 
just increasing the molecular mixing with no change in the mean profile would be a good 
start. The fact that the devices tested so far have not significant differences in the mean 
temperature profile does not necessarily mean that they did not produce any changes to 
the flow field. 



Figure 10 Figure 11 

Figures 7-10: Flow disturbance generators. 7: Serrated disturbance generator 1, attaches 
to inlet of the main pipe, 8: Vortex generator, attaches to inlet of the main pipe, 9: Swirl 
generator, inserts into inlet pipe near its exit, 10: Serrated disturbance generator 2, inserts 
into inlet pipe near exit, 1 1 : Spherical disturbance generator, inserts into flow path 
upstream of main pipe entrance 
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Figure 12: Temperature profile at the exit plane with and without the disturbance 
generator shown in figure 7 


Conclusions and suggestions for further work 

Efficiency is an important area of current hybrid combustion research. A cold flow 
experiment and a related numerical simulation motivated by hybrid rocket efficiency 
were described in the preceding paragraphs. The experimental and computational results 
do not match very well, and this issue remains to be resolved. An LES calculation or 
possibly a different turbulence model may be needed to improve the agreement between 
the two. 

There is a slight ambiguity in relating results from the cold flow experiment to the actual 
hybrid motor. While enhanced mixing is in general related to an improvement in the 
efficiency of a chemically reacting system, in the case of hybrid combustion, better 
mixing may result in a higher combustion temperature, which will increase the regression 
rate. If the cold flow experiment shows improved mixing for a particular device, it is hard 
to say if a similar device in the real motor would result in higher efficiency, or if it would 
instead increase fuel production at the fuel surface, in which case the efficiency could 
remain unchanged or even decrease. However, despite the ambiguity, it is worthwhile to 
do testing with the cold flow experiment. If a device does not produce significant changes 
to the flow field in the cold flow case, it is unlikely to affect the flow field in a real motor 
either, and it is probably not worth the effort and resources to test a similar concept in the 
combusting environment. On the other hand, if a device substantially alters the mixing 
behavior in the cold flow case, it is more likely to cause some type of a change to the 
flow in the motor. However, further testing in the reacting flow will be necessary to 
determine the precise effect. 



If it turns out that influencing the flow by upstream disturbances can not produce 
significant enough changes to the flow field, one option is to try to influence the mixing 
downstream of the port. If this idea is pursued, a section corresponding to a post- 
combustion chamber could be added to the cold- flow experiment downstream of the 
main pipe. In this case, there would be a more direct correspondence between the results 
of a cold flow experiment and the hybrid rocket efficiency. If the flow in the port is 
unchanged from one case to another, the fuel production at the wall should not be 
influenced. Enhanced mixing in the cold flow case should correspond to higher efficiency 
in the reacting case. 

Another possibility to be pursued in the future is addition of flow visualization to the 
experiment. Seeing the flow patterns may help in the design of new disturbance 
generating devices. 

Although the experiment has not yet produced positive results in terms of finding a way 
to mix the axial flow with the radially injected fluid more effectively, future efforts may 
turn out more fruitful. None of the devices tested have produced large changes in the 
mean profiles. However, some measure of molecular mixing should also be obtained for 
assessing their effectiveness. If future work is pursued with this experiment, it is 
important that this measurement is added. 
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Transient Modeling of Hybrid Rocket Low Frequency Instabilities 

M. Arif Kandicyogh*, Skne De Zihn 1 , Brin CntwdP ad GregZiffiac* 

Abstract 

A comprehensive dynamic model of a hybrid rocket has been developed in order to understand and 
predict the transient behavior including instabilities. A linearized version of the transient model predicted die 
low-frequency chamber pressure oscillations tha are commonly observed in hybrids. The source of the 
instabilities is based on a complex coupling of thermal transients in the solid fuel, wall heat transfer blocking 
due to fuel regres si on rate and the transients in the boundary layer that forms on the fuel surface. The 
oscillation frequencies predicted by the linearized theory are in very good agre e me n t with 43 motor test 
results obtained from the hybrid propulsion literature. The motor test results used in die c o mp a ris on cover a 
very wide spectrum of parameters including: 1) four separate research and development programs, 2) three 
different oxidizers (LOX, GOX, NjO), 3) a wide range of motor dimensions (Le. f rom 5 inch diameter to 72 
inch diameter) and operating conditions and 4) several fuel formulations. A simple universal scaling formula 
for the frequency of the primary oscillation mode is suggested. 


1) Nomenclature: 

k: 

Blowing parameter exponent 

A : 

Pre-exponential coefficient 

J : 

Laplace transform of oxidizer mass 
flax pertivbation 

Apy A m - 

Port and nozzle throat areas 

L>L m z 

Length of the fed port and motor 

B : 

Blowing parameter 


rj4. 

C: 

Specific heat of fed 

A, . 

c': 

Boundary fayerdefaytane 


Oxidizer and total mass flow rale 


coefficient 

a: 

Mass flux exponent 

* Q : 

Skin friction coefficient and Stanton 

O/Fz 

Oxidizer to fed ratio 


mmdxr 

p. 

Averse chamber pressure 

* * 



^ _. , . , , „ , 

» C dmo 1 

Measured and ealciibicd 

Rr- 




r 9 Rz 

fluxes 

Dimensional and non-dimensional 

D : 

Average port diameter 


E.- 

Activation energy 

Kc z z 

regression rate 
Reynolds number 


Eaogy hmefen 


Specific and univeisal gas constants 

/ * f\-L * 

Primary oscillation frequency, first 
acoustic mode frequency 

RT^z 

Average gas constant temperature 

F t ‘ 

Thermal system transfer function 

s z 

product 

Laplace transform variable 

Fjc- 

Thermal-Combustion coupled 

1z 

Time 


system transfer ftmrtiou 

Tz 

Tenycratiue 

G a ,G,i 

Oxidizer and total mass flaxes 

u,v: 

Axial and normal velocities 

K K- 

Total and latent heats of gasification 

U*z 

Diffusion speed 
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V V * 

P' r * ■ 

Port and motor volumes 

z z 

Axial distance along the port 

Sz 

Boundary layer thickness 


Chamber pressure oscillation 

» - 

ampmuoe 

7 ' 

Ratio of specific heals 

k z 

Thermal diffusivily of fuel 

Pf* P ~ 

Fud and average gas densities 


Boundary layer response time and 


thermal tag time 


2 : 

Gas phase response coefficients 

P- 

Average gas viscosity 

az 

Ambient 

I: 

First perturbation variable 

ez 

Free stream value 

ref : 

Reference quantities 

s : 

Superscripts: 

Surface 

- 

Non-dimeBsiaBal variable 

A- 

Laplace transfonn variable 


2)IMcte: 


The aikmaMe thrust oscil la tion s for an 
operational propriskm system at Smiled by ce r ta in 
practical considera tio n s induing a ccel er a tion 
loads on the vehicle stnactares and also on the 
payload. Has makes combustion sta b il ity a critical 
issue during the d evelopm ent of a near p rop u lsion 
system. The comb u s ti on instabilities of solid and 
liquid systems have been studied in depth trad 
partially s uccess ful theories along with c cr ta m 
practical rules for designing stable systems have 
been developed 1 ' 23 . Even with this extensive 
knowledge base, development of stable rocket 

resources during the design and testing phases of a 
new system. Even though he transient opera ti on 
and instability mechani s m s of hybrids have not 
been e xp lored as extensively as the more ma t ur e 
chemical systems, past experience shows that 
hybrids do not have the cat a strop hi c instabilities 
that liquid and solid systems may potentially 
present 

In fact, hybrid systems typically show 
finite amplitude (Le. 2-20% rms of mean chamb e r 
pressure) low-frequency cham ber pressure 
osaMons 45 ^’ 1 . This roost common hybrid 
instability is in the form of limit cycle oscillations 
with frequencies much ™IW than the first 


IcngrtndiBd acoustic mode of die combustion 
chamber. 


Even though the exact cause of die low 
frequency hybrid nwla hiH tirs is not yet known, a 
few piausMe t h eories exist 534 * 4142 - Most of the 
suggested theories lack die mathematical fo r malis m 
that exists in the solid and liquid instability models. 
The fundamental difficulty hi producing a 
mechanism that would generate instabilities comes 
from the fact that hybrid burn rates are typically 
independent of the rhambrr pressure. This, for 


the thermal lags in the solid and the chamb er gas 
dynamics which is die so urce of the L* i nstabi l i ti e s 
in solid rockets 1 . 


One particular theory, the T€ coupled 
model 1134 , which is based on die coupifrig of 
thermal tags in the solid and the lags in the gas 
phase transport, has predicted (based on a 
math e mati cal formalism) fuel mass production 
oscillations that would drive the o bserv ed p re ssure 
oscillations. The purpose of this paper is to extend 
the TC coupled theory to include the effects of the 
c h a m ber gas dynamics and co mp a re the theory 
predictions with the motor test data. The merhs and 
shortcomings of the TC coupled model will also be 
discussed. 


3)Tr 

Motors: 


i ■ Hybrid Racket 


Our ap proa ch to modeling the dynamic 
behavior of hybrids is to isolate the subsystems of 
the motor and consider every smgle subsystem 
mdrviduaily. For a full de scri p t ion of motor 
transi ents one has to consider the following 
subsystems. 

a) Feed System: In a hybrid motor, the liquid (or 
gaseous) oxidizer needs to be fed in to the 
combustion chamb er through a feed system. In 
reality the feed system response time is finhc due to 
the capacity of the eiemetes in toe system. Accurate 
response depends on foe details of the system, 
which is likely to be significantly d i ff e rent for 
every design. For this reason we will bypass the 
feed system dynamics in our mvestigations by 
assranmg the oxidizer mass flow rate as the input 
parameter. This is a reasonable assumption for most 
systems because isofating dements (Le. cavitating 
venture, sonic orifice) are present. 

b) Vaporization of the hkpM Oxidizer: 

Complete vaporization of the oxidizer droplets in 
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the combustion chamber requires a certain 
characteristic time depending cm parameters such 
as the droplet size and the thermal/flow 
e nvironm ent of the space where droplet 
vaporization takes place. For the purposes of this 
paper we will assume that Ihe oxidizer is in gaseous 
phase when it enters the port. 

c) DifladoB and Combustion in the Boundary 
Layer: It requires some tone for die hybrid 
boundary layer properties to adjust to the changes 
in the port velocity or the fuel blowing from die 
surface. The associated dynamics has been 
discussed in detail in Refs. 13, 14 and wiD be 
briefly covered in section 5 of this paper. 

d) Thermal Response in the Solid Grain: A 
change in the wall heal flux to the hybrid fuel grain 
can not be followed immediately by fuel production 
due to the futile heal capacity of the solid fuel. 

e) Chamber Gas D yn s nurr The c ha mber pressure 
responds to the changes in the mass flow with a 
time scale proportional to foe filling time of the 
cha m b er. Acoustic response is also resolved in dus 
sub-system. 

In foe snbseqnent sections foe last force of 
these sub-systems (sections c-e) wffl be m odeled 
and invcstijgaled individually. Eve n tual ly, these 
subsystems wiD be coupled to give foe overall 
system response. The order of magn i mrir estimates 
of foe tune scales of some of the important 
processes enrounlei cd ia a typical hybrid motor are 
lisled in Table 1. In our no dcl n g , solid and 
gaseous kinetic times are asisumrd to be fast 
compared to the other relevant time scales. 


4) Thermal Lag Model: 

The regression rate of the hybrid fuel grain 
cannot respond to the changes in foe surface heat 
flux instantaneously due to the finite thermal 
conductivity of foe solid fuel. A complete transient 
model for a hybrid rocket system requires a 
dynamic model to predict die regression rate history 
for a given wall heat flux schedule. In order to serve 
this purpose a thermal lag model that models the 
nonlinear heat conduction in the fuel gram with a 
moving boimdaiy has been formulated and reported 
in Refs. 13 and 14. The schematic of this particular 
model is shown in Figure 1. 

In die thermal lag model, die gasification 
and pyrolysis reactions at the surface are both 
modeled by an exponential expression of foe 
Arrhenius type. For chemical reactions, the 
exponential constant is am activation energy, 
whereas in vaporization it is the latent heat In order 
to describe this behavior, we assign an average 
effective activation energy, E a , resulting m: 

r-AeW- ( 1 ) 

Here T s is the surface te mperature , r is 
the regression rate and R* is foe universal gas 
constant. The mathematical formulation and various 
solution techniques for the thermal lag model are 
presented in Ref 13. One solution of interest is 
obtained by perturbing foe full nonlinear system 
around die nominal operating point The linear 
initial-boundary value problem defined for the first 
order perturbation quantities has been solved with 
use of foe Laplace transformation technique. The 
transfer function between the regression rate 
perturbation and foe applied beat flux pert ur ba ti on 
can be written as 


F *.w 

r aw t +V4*+l) )-2E^ +2E £ E^* 


( 2 ) 


Note that J^(x) and £,(,) are the 
Laplace transforms of the am-diroeusional 
regre ssi on rate perturbation and the wall heat flux 
perturbation, respectively. Here foe following non- 
dimensional parameters are defmed for 
convenience. 


E £ - AF 

E - ^ 


L CAT 

*r=(r,UT a 

Q= 9 W 

Qr*f 


Qr*f -r^PjCAT 
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to investigate the stability character of the TC 
coupled system. The co mm o n ly used method is to 
map the poles of the transfer function in the s ir 
plane. Particularly, the real component of a certain 
pole of a transfer taction indi c ates amplification 
rale associated with that pole. Similarly the 
ima g ina ry part rep r esents the oscillation freq u e nc y . 

First a system with zero boundary layer 
delays (Le. r M = r m = 0 is determined to be 
always stable. In fact plot of the transfer function 
for a system wife typical parameters taw s no 
poles in the s plane (only a zero at (0, 0)) indicating 
no sign of i ns tabi lity (Le. See Figure 3). If a 
positive delay is introduced be tw e e n the regression 
rale and wall heat flux, a scries of unstable poles is 

is shown in Figpe 4. Note that ail the other 
parameters are kept i dentical to the case with no 
delays. Even though their e xi s ts an Mimic scries 
of pedes g enera t ed (at die same amplif ica tion rale) 
we will only concentrate on the pole with fee 
lowest frequency (tadamental mode). We believe 
that the higher order inodes are an artifact of using 
a simple delay instead of the full dynamics for fee 

frequency modes, even if they exist in a real 
system, are fikefy to be tinu ped more effectively 
co mpared to the fh mlanw itil mode. 

ft is impor t «a to identify fee necessary 
conditions that must exist for fee ge neration of TC 
coupled instabilities. A careful eaanoamti ou of fee 
dcnannmOor of fee transfer inaction, Eq. 6, shows 
that for fee inst ab il it y to exist Mow i n g e xpo ne nt, 
k 9 activation energy, E a and delay time, 
must be nonzero. Hus indicates that fee instability 
is a result of the cosqdmg of three physical 
phenomena: Mocking of heat transfer by radial 
infection of fuel ma s s , t hermal tra n sients m fee 
solid and b o u n da ry layer dynamics. The gas 
phase/sotid phase coupfi r g m e chanism ca pt u red by 
the TC coupled model is shown in Figure 1. 

Next we consider fee effect of various 
parameters on fee TC coupled instability, ft can be 
shown feat fee effect of fee regression rate mass 
flux exponent rv, heat of vaporization of fee feel, 
L r , the beat capacity of fee solid fuel, C, and fee 
surface temperature, T s , on fee oscillation 
frequency and amp l i t u de are negligible for fee 
range of these parameters that are commonly found 
in hybrid systems. Also TC coupled system 
characteristics are completely independent of fee 


solid density, pj and fee heat diflusivity in fee 
solid, ir. 

The effect of fee activation energy on fee 
amplification and frequency is shown in Figure 5 
for T m =38 miltisrconds and k =0j68. As 
mdiraard by fee figure, for activation c uagi es 
co mm o nly observed m hybrids, ranging from 5 
kcal/mole to 60 kcal/mole, the variation in fee 
osciBation frequency is relatively small. The 
amplification increases with increasing activation 
en ergy. For tins specific case fee systems with 
activation energies la r g er than 3 kcal/mole show 
positive amplification and unstable behavior. 
Decreased stability at bigger activation energies is 
expected since at high E a values the surface 

te m per atu re and regression rate are more closely 
coupled. This coupling is one of the necessary 
ingredients for the g ener a ti on of the instabilities. In 
fact in the e xtr eme case of £ a =0, for which fee 
strface temp er ature is completely independent of 
the regression rate, oscillation frequency goes to 
infinity and fee amplification goes to negative 
infinity. 

The activation energ ies for fee polymeric 
systems feat are typically used in hybrids ranges 
from 10 kcal/mole to 60 kcal/mole. For fee 
Donpofymeric paraffin-based fuel system fee 
activMon energy is equivalent to the latent heat of 
v apor iz a ti on 1 *. Hus value for fee paraffin -based 
fuel formulations is estimated to be 17 kcal/mole 
(Ref. 13 )l 

Similar results can be obt ai ned for fee 
effect of the Mocking exponent, as shown in Figpne 
6 for r U2 = 038 and £*=15 kcal/mole. For fee 
range of values reported in fee literature for k (i.c. 
O-68-0.77) 15J *, the effect of Mowing exp onent on 
fee frequency is negligible. The amplification 
increases with increasing k and systems wife k 
values 0.45 and larger showed unstable behavior. 
Note that for fee physically unrealizable case of 
4=0 the instability disappears, since one of fee 
n e cessa r y coopting m e chanism s is eliminated. We 
finally like to note feat even though fee blocking 
effect in paraffin based systems is somewhat 
reduced due to the two phase character of fee flow 
field, its is still significant enough to estaMish fee 
coupling between the gas phase and solid phase. 

The most MlncntM parameter on the TC 
coupled frequency is determined to be the 
boundary delay time, r^. Figure 7 shows the 
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predicted frequency as function of the delay time 
for three activation energy values, E a = 5, 15 and 
50 kcal/mole. As shown in the figure, the effect of 
the activation energy on the frequency, especially 
aft longer delay tones, is smriL For Ae purposes of 
this paper the effects of activation energy and 
Wowing exponent wiS be ignond and Ac 
frequency will be represented as a function of the 
boundary delay time done. The JoDcwing crave fit 
equation is to quantify the inverse 

rela ti on shi p b e t w een the frequency and Ac 


/- 


0.48 

T tn 


( 7 ) 


As shown in Figme 7, this e xp re ss ioB 
represents the p re di cted frequency to a very high 
degree of accu ra cy and will be used for all Ae 
hybrid systems that will be discussed later in Ais 
palter. 

7) Gas Ifyawir Madrtr 


In Ae previous section the d ynamic 
models for Ae thermal lags and Ae hybrid 
co mbust ion are summarized. In addition, these 


subsystem models are coupled 
respon s e of the fad wmk 
in the mass flux m the port. Howe v e r, ] 
significant parameters of Ae jw mhrt operation such 
as the c ha mber p ress ur e, sp e cific i mp ul se and 
thrust cannot be obtained solely from the TC 
coupled system. These variab les can only be 
determined after introducing a model for the gas 
dynamics which will be discussed in Ais section. A 
gas dynamic model, in general, should use the 


oxidizer mass flow rale and the fuel mass 
ge nera ti on rate as therapists and it should yield the 
param e te rs that are more closely related to Ae 
p er fo rmance of the motor such as Ae cha m be r 
pressure and Ae motor O/F ratio as Ae outputs^ 


fat onr specific model (2V-port model) we 
divide Ae chamber into three on pcaafts, pre- 
combustion c ha mber , port volume and the post- 
combustion cham ber (Le. see Figure 8). We treat 
the pre and post com bu s ti on « lumbers as zero 
d im en sio nal volume ele m en t s , but consider Ae 
variation of stale variables along Ae axis of the 
port. The mathemati cal fbnnulmion is obtained 
after the application of the conservation laws for 
the three volume ele m e n t s. The set of equations 
resembling our model is solved analytically after 
l in eariz ati on and as a result of that operation, a 


transfer function for Ae gas dynamic subsystem is 
obtained. Some numerical simulations are also 
per fo rmed to determine the validity of Ae 
approximations used in Ae perturbation solutions. 
Note that a detailed description of Ae gas dynamic 
model is given in Ae Appendix of this paper. 

The following co n clusio ns can be down 
from Ae resufts of Ae gas dynamic model as given 
in detail in Ae Appendix of Ais paper and in more 
detail in Ref. 13. 


Gas dynamic model ca p tures the 
filfaq^emptying dynamics and also the 
longitndnial acoustic respons e of Ae hybrid 
combustion chamber. As can be seen from 
Figure 13 a cousti c modes are stride (well 
damped). No unstab le oscillation modes are 
generated by Ae ^s dynamic system by hsdf. 
The nonlmear ifishibances exche Ae chamber 
acoustic modes. Snnriations showed that a 
sudden pulse in Ae oxi di z e r mass flow me 
excites Ae a c o us ti c s modes. These anodes 
decay in a relatively short period. This ringing 
p henomenon shows that Ae longitudinal 
acoustic modes can be driven by nonlmear low 
frequency oscflbtkns. Namely the nonlinear 
low fr e qu ency p re ss me waveforms 
continuously excite Ae acoustic modes of Ae 
cham ber . AMROC motor DM411 showed aft 
least 5 longitudinal acoustics modes. All those 
were wefl dam pe d. The gas dynamic model 
when c orrect e d for isentropic speed of sound 
predicts the acoustics frequencies of the DM- 
01 motor to a high level of accmacy (See 
Table 7). 

Gas dynamic model coupled with the 
boundary layer dynamics (r M * 0) does not 
produce any ins t ab ilities (no extra poles are 
generated in Ae s plane). This conclusion is 
based on transfer function analysis and also 
mmcrical sreratokras. This mqmrtst result 
highlights Ae importance of the thermal fag 
system (ix. and its unique dynamic capability 
of p rodu ci n g a phase lead) in p rodu cing a 
co up li ng me c hanism to generate low 
frequency unstable modes. 



In this section we integrate Ae gas 
dynamic system to the TC coupled system to 
establish Ae TCG coupled response (Le. See 
Figure 2). The TCG coupled system yields the 
performance of the rocket (such as diamber 
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pressure) for some given input of oxidizer mass 
flow rate, winds is a fundamental control 
parameter in hybrid rockets. It is fair to state that 
the TCG co upled system represents the most 
fundamental dynamic behavior of a hybrid rocket. 

The s ch e ma tic of the hybrid subsystems 
and the infor mat ion flow between the subsystems 
in the context of TCG coupfing are shown in 
Figure 2. Note that the lower block in the 
schematic re p rescids the TC coupled system. The 
input for the TC coupled system is the local mass 
fhxx info rma ti on and the o u tp ut is the mass 
generation or Ac regression rate. The upper block 
shows the 2V-porl gas d ynamic modcL The gas 
dynamic model takes die oxidizer flow rate and the 
mass genera ti on rate information and reveals the 
imp or tant p er f orm an ce p a ram eters such as the 
chamber pressure and specific impulse. The TCG 
coupled system can be considered as an overall 
model for the d yna mi cs of a hybrid propulsion 
system that utilizes gaseous coocfizer wfaicb is 
delivered by an isolated feed system. 

The generalization of the TCG coupled 
system to a liquid hybrid with si gni fica nt feed 
system d ynamics can be ach ieved easily. Han dy 
two more modules for the fitprid droplet 
evaporation and the feed system dynanics must be 
added to the front end of the TCG model. We 
believe that die fundamental phenomenon 
ge nera tin g the low frequency instabilities of hybrid 
systems does not involve either of the evaporation 
lags or die feed system dynamics. 

In the TCG coupled system the gas 
dynamic module converts the fuel mass oscillations 
produced by the TC coupled system ido the 
chamber pressure oscillations, h has been 
determined that m the p ro ce s s of conversion, the 
oscillations p rod u ced by the TC coup l ed system 
(ix. frequencies and nullificati on rales) are not 
altered. This fact is de mo ns tra ted in Figure 9, 
winch shows a plot of die TCG system transfer 
function for the TC coupled system depicted m 
Figure 4. The same plot also indicates poles 
associated with die hflmg/emptying mode (find 
order nou-oscdlatoiy mode) and dm longitudinal 
acoustic modes with negative read co m p o n e n ts 
(damped). TCG coupled system captures all the 
necessary mderiying dynamics ( lincm iz ed verson) 
of a gaseous hybrid with a decoupled feed system. 


9) NASA Ames Paraffin-Based Motor Tests: 

A promising fast burning paraffin-based 

fuel has been tested extensively m the Hybrid 
Combustion Facility (HCF) at NASA Ames 2 *. The 
motor outside diameter is 10 inches and gaseous 
oxygen is the oxidizer for all the tests conducted at 
the HCF. The p re ssure time trace of a typical 
paraffin-based hybrid motor test is shown in Figure 
10, whereas the Fourier transform for fee same test 
is given in Figure 11. Note that the chamber 
p re s s u re traces for classical hybrids as presented in 
the lite rat ure are very similar in nature. As 
indicated in F%nre 10 for a particular case, die feed 
system and fee cha mber were completely 
decoupled for a! paraffin-based motor tests used in 
this paper. 

The Fourier transform shows three broad 
peaks conesponifing to fee low frequency hybrid 
instability, fee Helmh oltz mode and fee first 
longitudinal mode. It has been generally observed 
feat fee dominant mode, low frequency oscillations 
are a cco mp ani ed by lower amplitude higher 
frequency acoustic modes. See Ref 6 for an 
extensive discussion of fee stability c h a r acte ri stics 
of fee paraffin-based motors. 

Motor Test Data 

In this section we wifl compare fee TC 
coupled theory osefltation frequency predictions 
wife fee c hamb er p re ssu re oscillation frequencies 
from by brid motor test data. hi order to produce a 
conclusive c omparison, we have selected a large 
number of tests form several hybrid rocket 
development programs. The key features of each 
progr am is discussed briefly in fee following 
paragraphs. 

- AMROC Molar Tests 443 : 1b this case the 
reported resuhs arc for 4 di ffe rent thrust class 
hybrid motors which aB utilized LOXAJOX as 
fee oxidiz er and HTFB as the fuel. The data 
from fee AMROC tests is given in Table 2. 
HPDP 11 Inch Motor Tests 51 : The oxidizer 
used in these tests is LOX and GOX and fee 
fuel is HTTra/Escorez formulation. The motor 
case diameter was 11 inches and fee data from 
5 tests used in this study is given in Table 3. 
Only the tests feat exhibit comb u s ti on 
instability and p e rf ormed wife single port 
motors are considered. The runs wife fee 
rearward-faring steps at fee port entrance are 
also excluded m this study, since in these tests 
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the flow field in die port is significantly altered 
winch may cause a major cha nge in the 
bo unda ry layer delay tme co effici e nt . We abo 
include test #1 of 11 inch LOX motor data in 
die analysis. We only consider test 1 since it is 
die only ran with a c h am b er configuration 
identical to GOX motors used m the 

- JDtAD Motor Tests 1 ”': The oxidizer used m 
these tests is GOX and the fbd is 

HTPB/Escorez formulation. The motor case 
foameter was II inches and the data from 8 
tests used m this study is given in Table 4. 
Only the tests that exhibit combustion 
instability and p er fo rmed with smgle port 
motors are considered here: 

Arizona State University 22 : One test is 

included in the com pari so n for a 5 inch 

d ia me te r nitrous oxide motor. The fuel was 
HTPB (Le. See Table 5). 

Paraffin-Based Motor Tests 2 *: These tests are 
conducted at NASA Ames in a 10 inch 
diameter test facility. The oxidizer was GOX 
and the fnd is two separate nonpoly menc 
paraffin-based formulations (SP-la and SP-4)l 
The data form 25 motor tests are r epor te d in 
Table 6. 

The prediction of foe frequ en cy requires 
the estimation of foehoondary layer delay time for 
each test from foe reported data. Since foe reported 
data is some wha t differen t for each p r o g ram 
d iffe ren t scaling laws for foe delay time must be 
derived. 


The first one of those is for foe various 
size AMROC motors which operated at different 
L levels. Since the information on the AMROC 
motors that can be found in the literature is limited 

to L , c and motor O/F ratio, it is desirable to 
express foe boundary layer delay time. Equation 4 
(with z- L/2) in terms of those variables. The 
average velocity in foe port can be approximated as 

, GXl+lO/FW+O/FpT^ 

mr 2P C ' 

where RT m is an average value in foe port. 
Substituting this expression in foe delay formula 
aid using foe relations for foe total mass flow rale 
m t = G,A p mid port volume V p — LA p yields 


With the use of the total mass flow relation, 
m t — P c A m j and foe definition £ = V m fA u , 
the delay equation can be written as 


f bi2 


r, ft+O/FW+M/Fg , . 

V m RT" 


Here V p jV m is foe ratio of foe port 

volume to foe motor voftsme, which is estimated to 
be appnreinalety 0.8 for AMROC DM-01 motor. 
We assum e fort afl AMROC m o to r s posses s foe 
same average gas co n sta n t temperature product, 
RT^ aid volume ratio, V p jV m . After noting fort 

all motors operate at very sanilar O/F ratios and 
values, it can be stated fort in foe series of 
tests the boundary layer characteristic delay time is 
proportional to I*. 

In order to reduce data form the rest of foe 
tests, a sanilar relation for foe botmdsny layer delay 
tone in terms of foe cham ber p re ssu re and port flux 
levels must be derived. The port velocity can be 
replaced by an average value 
with given by foe gas law P c = p mr RT mr 
to yield finally for the hown fa ry layer fag: 






(G„ + G,Ver„ 


<ii) 


Equation 11 gives the scaling law for foe 
boundary layer delay in terms of foe oper atin g 
points of foe rocket motor and also with foe size of 
the motor. This states that delay increases with 
increasing c ha mber pressure, port length and 
decreases with increasing mass flux in the port. An 
im p ortant conclusion that can be drawn from those 
observations is fort an i n cre a se in the dumber 
pressure or a decrease in foe port mass flux act to 
decrease the oscillation frequency of this type of 
TC coupled system. 

The total mass flux can be efawwnatrd form 
Equation 11 by introducing the average oxidizer to 
fuel ratio of the motor. 



l U2 


V _Zl 

m,Klf20/FV(l + 0/F)Jrr„ 


(9) Equation 12 will be used to reduce the paraffin- 
based motor data. 
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Note that equations 10, 11 and 12 are 
equivalent and the constant muhipiier c is 
identical for all of these expressions. 

The average temperature gas constant 
product that a pp ea r s m the denominator of al the 
boundary layer defay fame eq u ations is assumed to 
be constant for til of the LOX/GOX motor tests 
that is con sidered m lias paper. This is a fairly 
good as s ump t i on since RT is a weak fimcfaon of 
the motor Q/F for the practical opera ti ng 
conditions 13 . In fad the m aximum expected 
variation on the RT values for all the LOX/GOX 
tests used for comparison is predicted to be less 
than 5%. The a b solute value of RT {rwc is selected 
to be 639 10 5 (m/sec) 2 for LOX/GOX motor tests 
and since the same value is used for aO calculations 
any error ia RT^ wifl only effect the munericti 
value of the em p i ri c al defay cons tan t c . 

In order to double check the RT^ value 
used in the calculations one could use the fast 
acoustic frequency m ea sme d form the motor test 
data and use dial to estimate the speed of souad 
and finally RT^ .based on the following equation. 


Here y is the ratio of the specific heats averaged 
over the length of die motor, f^_ L is the first 
longjtndinti acoustic mode and is the effective 
length of the motor. For a more accu r ate 
calculation the c om p lex gas dynamic model 
discussed in die Appendix of this paper can be 
used. AH of the NASA Ames motor tests with 45 

arorrad 370 Hz. For these tests with use of y— 1.25 
and 1^=127 m, one obtains RT^ of 7.1 I0 5 
(m/sec) 2 which is only 10% higher than the 
assented value. RT^ for the NjO test is selected 
to be 5.1 10 s (rn/scc) 2 . This is 20% lower than die 
RT^ for the GOX/LOX system due to the low 
t e mper a t u re of die NjOAydnxaiboa c o m b u s ti on 
p roducts. 

Now Equation 7 can be coupled with the 
expressions for the boundary delays time to 
evaluate the frequencies predicted by the TC 
coupled model. 




Jr 

2L m 


(13) 


/ = 0.48- 


RT_ 


’{v p /V m ) K|+O/fV0 + 2O/^ Lc 


(14a) 


/ = 048 ^ G ° +G ‘ ^ 7 ‘"’ r = 0.48 
c'LPr 


f 2 » K^= 

( O/F) cLP t 


(14b) 


The oscillation frequency as a func ti o n of 
the bomdary delay time predicted by the TC 
coupled theory along with data form the hybrid 
motor d evelo p ment p rog ram are p lotte d in Figure 
12. Note that we have used a bowdary layer delay 
constant value, c\ of 201 for til motor tests m 
order to obtain the best fit be t w ee n the theory and 
test results. Please note dial this e xp erimental vadac 
for c is the same order of ma g ni t ude with the 
rough estimate (Le. -035) for a bounda r y layer 
with no blowing and do combustion. 

Considering the wide range of test 
conditions, the agreement between the measured 
pressure oscillation frequencies and the TC 


coupled model predictions is surp r isin gly good. 
Note that the comparison inefades 43 motor tests 

• Using resuis from 4 diffe r ent 
development and research programs 

• Using 3 separate oxidizers (Le. LOX, 
GOX and NtO) 

• Covering a wide range of motor 
dim en si on s (Lc. from 5^00 Arizona stale 
motor up to 72" OD AMROC motor) 

• Uring many different fuel compositions 
(Le. HTPB, HTFB/ESCOREZ, 
HTPB/DDI/ESCQREZ, 2 paraffin-based 
formulations) 
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• Covering a wide range of operating 

conditions. 

Equation 14a radicates that for motors 
operating at the similar O/F, the oscillations 
frequency is inversely proportional to the L of die 
motor. This conclusion is consistent with die 
previous results reported in the htendrae. Similarly 
Equation 14b shows for motor operating at equal 
oxidizer mass fluxes, chamber pressaes and O/F 
the oscillation frequency decreases with increasing 
fuel length. In order to pro ve that motor 

length plays a critical role setting the oscillation 
frequency, we have plotted the Ames motor data, 
which contains two sets of grain lengths, by 
removing the length effect horn the seating law. 
This resulted m an increase ra the scatter in the 
frequency data mdka ti ag die i mpo rta n ce of the 
grain length A similar check to co nfi r m the 
im por ta nce of length can be conducted by 
comparing die JIRAD and HPDP motor tests (Le. 
ran with g ain s 102 inches long) with the Ames 
p ar affi n tests (Le. ran with grams 32 or 45 inches 
long) that were conducted under very similar 
opera ti ng conditions^ For example co m pa ri ng test 8 
of the HPDP program to paraffin based motor test 
4P-01 shows hat gram length must be included in 
the scaling law m order to explain the discrepancy 
in the o scilfali on frequencies for th e se tests. The 
same conclusion can be drawn by comp arin g test 2 
of HFDP p rog ram to the test Ames test 4L-08. 

TC coupled theory predicts oscillations of 
the r egression rate/fud mass generation rate and 
would he observed directly in a hybrid motor as the 
radial oscillations of the diffusion flame within the 
boundary layer. The gas dynamic model converts 
the fuel mass g e ner a tion rate oscillations to the 
chamber pressure oscillations, h is expected that a 
system operating at higher O/F ratios would 
produce less c ha mber p re ssu re oscillations since 
fuel mass constitote a smaller fraction of the total 
mass expelled from the nozzle. This fact is 
observed in NASA Ames motor tests, namely high 
O/F motors, in general, presented lower amplitudes 
c hamb er p re ss ure oscObtkras compa re d to the 
mean cluunher pressure. 

The simple linear model presented in this 
paper successfully explains the gener a t i on of small 
amplitude low f requency oscillations and the 
scaling of Ac frequency. But it also falsely predicts 
an indefinite growth of chamber presage 
oscillations, bn reality, the nonlinear mechanisms 
that exist in a hybrid motor would timh the 
indefinite growth of the oscillations and result in 


limit cycle oscillations that are commonly observed 
in motor tests. Moreover the simple theory does not 
explain when these low frequency instabilities will 
lake place and what then amplitudes would be. h 
has been shown in the hybrid development 
p ro grams that the fore end configuratioti/votame of 
the motor phys a paramount impor ta nc e in setting 
the amplitude of die low frequency inst a bi lit ie s . 
The motor tests indicated that, at least for gaseous 
motors, axial injectors result ra more stable 
operation co m p a red to the radial or conical 
mjection of the oxidizer. 

A natural e xplanati on for these 
observations is drat the TC coupled instability 
me chanism almost always exist ra hybrid systems, 
bat the amplitude of the limit cycle oscflbtions will 
depend on the excitation level by the disturbances 
in die right frequency range. We believe that the 
fore end fluid dynamics (Le. vortex shedding) is 
the priuKuy sourer of disturbances in a hybrid 
rocket For exanqde cha n g in g the injection scheme 
of the oxidizer from axial to radial, may introduce 
a flow disturbance c om po ne nt hat is p r e fe rr ed by 
the TC coupled system. If one assumes that all the 
kvw frequency oscillations sue developed by the 
fluid mechanics a the fore end of the motor, one 
would fail to explain die length effect on the 
oscillation frequency. Also all oxidizer injection 

seem to produce the same oscillation frequency. 
This feet is also difficult to explain by the 
assumptions that all die oscillations are produced 
from the fluid dynamic everts at the fore end of the 
motor, since the injection scheme should alter die 
fluid mechanics sig nifi can t ly (Le. vortex sheddrag 
frequency). 

We would like to note that most of 
musical instru ments also work on the same 
principle 2 ’. For example in a floe organ pipe the 
column of air inside the pipe is set into v ib ra ti on 
(Le. at its natural frequencies) by an edge tone 
which is produced by a jet of ak feat impinges on a 
carefully designed frp. The vortex shedding at the 
lip sets the edge tone which itself excites the 
natural inodes in the pipe. Note that in the case of 
an organ pipe the jet speed and die lip geometry is 
carefully selected for die pipe length such that the 
edge tone produced at the lip cootaras a large 
component at the natural frequency of the pipe that 
needs to be excited. In this organ pipe analog, the 
lip is the fore end of the motor (Le. geometry, 
injection scheme, velocity etc) and the o r gan pipe 
is die hybrid system modeled by the TC coupled 
theory. 
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11) Conclusions: 


The following conclusions can be drawn 
from the results of the TCG coupled model: 


A plausible m r dwm s in that generates low- 
frequency chamber pressure oscillations is 
suggested. This mednmsm is <hie to the 
coupling of die thermal lags, the gas phase 
comb u s ti on and gas dynamic sub-systems of 
die hybrid rocket The physical transient 
model of the motor is presented in a 
mathematical formalism that allows one to 
wiintug the expected oscillation frequency. 
The amplitude of die oscillations can not be 
determined by this linear model. 

* The oscillation frequency estimated by the 
model is in very good ^ re e ment with the 
motor test data from several programs utilizing 
three distinct oxidizers, several fuel 
formulations, a wide range of motor 

sugges t the following equation as a universal 
sealing law for foe primary hybrid oscillation 
frequency: 


f =0.239 


f 2+ 4_]as«. 

I O/F ) L P c 


(15) 


RT^= 639 10 s (m/sec) 2 for GOXAjOX systems 
ST aK = 5.11 10 s (m/sec) 2 for low energy oxidizers 
systems such as N^O 

This mode of i ns ta bi li ty is common to all kind 
of hybrids: liquefying or co nventio nal fuels, 
liquid o xi d iz er or gaseous oxidizers, h is 
reasonable to believe that this low frequency 
mode is present in every hybrid system to 
some extent. Some motors are more unstable 
com pare d to foe others because these motors 
do possess more dis tur b an ces to excite foe TC 
cou pl ed oscillations. For example foe oxidizer 
injector config urati on or pre-combustion 
cha mber geo m e try affects the scale and the 
frequency of foe distur b an ces that would 
excite the commonly observed low frequency 
mode. For tins reason the design of the fore 
end of the motor has been critical in 
controlling the stability of foe system. 

Even though the linea r theory very 
successf ul ly predicts the os eflfafti on frequency 
of the cham ber pressure, it foils to establish an 
estimation process for the amplitude of the 
oscillations. The theory predicts an imfimked 


growth of the oscillations. In reality the 
amplitude of the chandler pressure fluctuations 
will be limited by nonlinear effects that are not 
covered in the linearized TC coupled theory. 

It has been shown that the longitudinal 
acoustic oscillations that are commonly 
coexist with the low frequency oscillations can 
be generated by the excitation induced by the 
nonlinear waveforms of the low frequency 
chamber disturbances, hi short, it is plaus&le 
that high frequency acoustic modes are driven 
by the low frequency oscillations. 
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14) Appendix 

In our gas dynamic model the motor is 
divided into four subsections: 1) pre-combustion 
chamber, 2) port volume, 3) post-combustion 
chandler, 4) nozzle. Each component is 
approximated a ppropriately with one or zero 
spatial d im en s ional models. The subscripts 1 and 2 
are used for the physical properties in the pre and 
post combustion chambers resp ec ti vely. Various 
stations of significance are also distinguished by 
use of subscript notation. Namely the entrance 
station of the gaseous oxidizer into the system is 
re pr es e n ted by a subscript o, the nozzle entrance is 
denoted by n and the port entrance and exit are 
shown by a and h, respectively- Our models for 
these subsections are discussed in the next section. 

Pre-Combustion Chamber 

We use a zero dimensional gas dynamic 
model for the pre-combustion chamber. Thus the 
momentum equation and die energy equations 
reduce to the mean stale in forma t ion of the 
pressure and temperature in this volume. But, the 
mass balance equation is nontrivial ami it can be 
written as 



Note that the size of the vohnne is 
assumed to be c obsM over the bun time. We 
farther assume a co n stant molar mass and a 

polytropic process, P cc p” r in the chamber. Here 
n p is the polytropic exponent We believe that 

these are reasonable assumptions for this relatively 
cold volume which a c co mmodates negligible 
combustion. Under these staled conditions and with 
the use of the idea! gas law, equation A.1 can be 
reduced to 

tiP V 

— L = m 0 -m a where a, = — (A2) 
dt npRyi^ 
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Port Chamber 

In our model we treat the port as a quasi- 
ID tube with cont inuou s mass and heat addition. 
The conservation of mass can be written as the 
following diffe ren tia l equation- 


dp d G 

vm f 

dt dz * 


(A 3) 


Note that for the simplicity of notation we treat the 
gas cons tant t e mpera ture product as a single 
dependent variable. 

Furthermore we assume that the gas 
constant tempe ra ture product at the exit of the port 
volume is a fraction of the nozzle entrance value 
RT C which can be calcubled from the equilibrium 
chemistry at the operating CVF ratio of the motor. 


Here is the mu generate d at the axial 

location z of the port. At tins stage we prefer to 
express the mass generation term in the general 
functional form. 

*»/ =/(G,z) 


KT 2 =f b KT c (A.7) 

In short, we specify the whole temper a t ure 
field ia tire motor rhamhrr for a given operating 
point and we i gnore the variations in the 
temperature field dini n g the tra n si e nts In that 
re spect 2V-port is an isot h e rmal gas dynamic 
moon. 


The momen tu m equation in the port can be written 
in toms of the local mass flux, G, the local 
pressure, P and tiie local density, p. 


dG 

~dt 



| 3P 

r CfC y 

G 2 

d z 

d z 

2^7 

P 


(A.4) 


i n cl uded to capture ihe effect of tire skin friction on 
the force balance. Here C f s tan d s for he skm 

friction coe ffici en t at location r, A p as die port 

area aid the C for the circumf e rence of the port 
All the geometrical p roper ti es are assumed to be 
independent of time and axial dimension This is a 
reasonable approximation for typical hybrids which 
posses a uniform port diameter distribution along 
the length and also a slow varia ti on of the port 
geometry 15 . 


Post-Combustion Chamber 

Similar to the pre^combustion chamber 
treatment, we as si gn uniform thermod y namic 
p rope rti e s for the gas in the post-combustion 
cham ber volume. Th erefore the dynamics can be 
represented with the mass balance which takes the 
following form after the fixed volume assumption. 



Here m k is the mass flow rate entering the hot 
volume and m m is the mass flow rate exiting the 
volume through the nozzle. Note that we have 
inhenntly assumed that post-combustion chamber 
behaves like a well stirred reactor with the 
reactions fast compared to the transient rates of 
concern. Equation A.8 can be written as 


We further assume tint the ideal gas law 
holds locally m the port. 

P = pRT (AS) 

The closure of the system requires the 
energy equation. For the sake of simplicity we 
replace the energy equation with the following 
linear variation of tempe ratm e gas constant product 
in the port. 

+((/fr), (A6) 


a 2 ^L=rn k -m m where a 2 =^r (A.9) 

We have approxima ted the gas co n sta nt 
t em per a t u re p rod u ct m the post-combustion 
cha mber with product at the port exit. 

Nozzle 

We first assume that the nozzle flow is 
quasi-steady which is only valid at very slow 
transients. We also assume that through out the 
operation of the motor, the nozzle is always 
choked, namely at any time the ratio of the 
chamber pressure to the ambient pressure is larger 
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than the critical value associated with the 
p ropert i es of the expending gas. After these 
sHnpfifyuig ass um p ti o ns the nozzle flow rate can 
be expressed with use of the classical quasi- ID gas 
dynamic formula 



Here j c is the ratio of the speci fic heats of the 
combustion products. 

At moderate frequency oscillations (such 
as the acoustic modes of the motor), the acoustic 
a dm itt a nce of die nozzle must be considered and a 
more sophisticated behavior must be mode l ed. For 
our purposes, it is adequsne to adapt Crocco’s 
expression 24 for the nozzle transfer function which 
is valid for small deviations from the mean 
c ham ber p re ssure and for moderate frequencies 


^=-= A ^J f * +k m — (A.11) 

P2 V»T * 

Here k m represents the effect of mertia of the gas 
m the convergent portion of the nozzle in causing a 
phase diffe ren ce be t ween die oscillation of the 
velocity aid p res su re at the entrance of die nozzle. 



This integration is p erformed on the 
whole convergent portion of the nozzle. Note that 
Cm is the critical sound velocity, v is die local gas 
velocity in the nozzle, » r is the gas velocity at the 
nozzle e ntr anc e. 

Equations A.9 and A.11 and the ideal gas 
law can be combined to relate the pressure 
perturbation directly to the mass flow rate 
p ert ur ba ti on at the exit of the port. 



(A. 13) 


Note that a fraction of the nozzle volume 
is added to the post-combustion chamber volume to 
account for the inertial effects of die gas in the 
nozzle. For typical hybrids the effective nozzle 
volume is small co mp are d to the post-combustion 
chamber volume. We Anally indicate dot for our 
isothermal treatment the tempe ratu re derivative 
term is identically zero. 

Equation A.13 includes some properties 
of the c ombus ti on products which are not known 
offhand. These p rope rti e s are calculated from 
chemical equilibrium consideirions with use of 
the code STANJ AN for the selected propellants. A 
sample calculation result of the product, RTc with 
the motor O/F ratio for HTPB-ootygjen system is 
given in Ref. 13. 

S oluti o n Techniques 

In this section we develop small 
perturbation solutions and also numerical solutions 


for the system of equations derived for the 2 V-port 
model . 

Perturbation Solutions 


We determine pertmbation solutions 
around a nominal opera tin g point of die rocket for 
the gas dynamic system of equations described m 
the previous paragra p hs. These p ert u rba t ion 
solutions will be used to obtain a transfer function 
rep r esenti ng the hear response of the gas dynamic 
system. The transfer function uses die oxidizer 
mass flow rate as die inpu t and the system pressure 
for the output The details of the derivation are 
shown hi Ref. 13. 

The transfer function specifically between 
the oxidizer mass flux hput and the pre- 
combustion chanfoer pressure outjw! can be written 
as 
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*?■ {* + K^7t h ^ / 

* ^+^ 2 ) * +e ] A+ l^ter "*)+2] *+]?LJ>+*-{ F n:- a +% (l^^+e^nl^AL/l) <A ' 


This function also represents the overall nature of 
the response for gas dynamic system. Similar 
expressions for the other variables can be also 
derived, hot tor the psrpose of our arguments this 
expression is adequate 

Numerical Simmiahom 

In order to check the validity of the 
assumptions used in the p er tur b a tion solutions, we 
also perform mmerical simulations on die 2V-port 
model described in the previous section. We start 
with the finite di fferc achg of the special derivative 
terms and use the second order central differencin g 
over a uniform mesh along the axis of the port 25 . 
For the time m a r chin g we use the Rnngp-Kufla 
fourth order method. The boundary conditions for 
the port equations are obtained from the pre and 
post c ombusti on diwtur dynamic equations. The 
numerical method used in the r a i m! a t io n s gave 
stable and accurate remits for a typical time step of 
3x1 OT* seco nd s and 80 mesh poiats in the port. The 
details of the nu mer ica l scheme is discussed in Ref 
13. 

The simulation code is applied to various 
transient events such as the throttling and the 
sinusoidal oxidizer flow rate case. Some of the 
results are discussed in the following section. 

The response of the gas dynamic system 
for a hybrid motor can be deduced from the gas 
dynamic system transfer faction. For example, the 
stability features of the system are revealed by the 

which is the characteristic e qu ation. However any 
investigation of the gas dynamic system response 


re q uires an inpu t for the comb u s ti on response that 
is represented by the F w term. 

In this section we use the sim p le s t quasi- 
steady combustion model. Namely, all the tra n sien t 
effects including die thermal lags in the solid and 
also the boonefany layer delays except the gas 
dynamics are ignored in this model. For that 
reason, this model which does not allow any 
dynamic coupimg between the sub-systems of the 
rocket is highly desirable to isolate the 
investigations to purely gas dynamic p h eno m ena. 
The combustion respon s e func ti on can easily be 
derived from the genera) regression rate law after 
linearization and Laplace transformation. 

F n = a t nG *•' (A. 15) 

For the case of die quasi-steady 
combustion the gas dynamic system transfer 
function is plotted it Figure 13. The operational 
parameters of this motor are discussed in Ref 13. 
The plot shows a pole on the negative real axis and 
a row of unagmary poles also in the negative read 
half plane. The real pole resem bl e s the volumetric 
cap ac ity of the cham ber, whereas the imaginary 
poles represent the fongitudhiat acoustic nature of 
the cavity. Note that all the poles have negative 
real parts, namely the both acoustic and also the 
global response ccmponats decay in time. The 
conclusion is that the pure gas dynamic system is 
stable m nature. 

In order to test the validity of the t ransfe r 
function we consider die trivial case of no mass 
addition and no temperat u re or velocity gradients 
in the port. In tins case the transfer function 
simplifies to 


e, 


fim 


“* M jL 1 +L 2 ) ,+*]+ 


RTJ 

2s 

for 


tanh| ter“"wW2] 


(A. 16) 
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After the application of the L’HosphaPs rule one 
obtains the wefl known transfer function for a 
simple single volume system. 


e,= 


RT C I 

(^i +^2 -r+e 


TfS + l 


(A.17) 


H e re Tj = L fe * is the charac te ris tic filling time 
of the motor chamber. The other parameters are 
defined as i* = VfA p and c =T c Jin^ whidi 

are commonly used in rocket ti te atoe. Note tint 
the total vohmeof the motor comb u s ti on cha mb er 
can he written as V = + L 2 + L)A p . This is the 

transfer function generally used (especially in the 
industry) in hybrid dynamic modeimg. 


The preceding equation does not predict 
the filHng/emptymg time scale of a hybrid rocket 
accurately. For than reason we develop a better 
estimate by the fo ll o wi n g argument, hi the region 
of the s plane that the pole for the fiflingfcxnptying 
process is expected to be located, the following 
approximations can be made. 


b * s0 ’ < A - ,8 > 

Under these approximations the gas dynamic 
tr an sfe r function simpl ifies to 


{l, +L 2 +| L, L 2 ±2/{F n - 5+e J- <A ,9) 

Note that since F^-a # lenab-dtpe-tatt of s, this transfer function represent the behavior of a first 
order system widi a characteristic time of 



z + ^ir“^l tanh^F^ -a* 

I 1 / 2 ) 

e 

jl-tanfaJ/^-o, 


1 


(A.20) 


The throttling re sp on se for a throttling 
ratio of 12 of the motor with the sam e 
specifications is calculated with d iifc i c nl methods 
and (dotted in Figure 14. Figure shows the time 
variation of the non-dimensional pre-co mbu stion 
chamber pressure defined as 

H'VrtW/ri) where Pj is At initial 

p re ssu re and is the final pressure. The striking 

observation is that the pert u rbat i on solution wftfa 
the fiDing time defined by equation AJ20 » in 
excellent agreement with the numerical simulation 
resulL This shows that equation A20 gives a good 
estimate for die filling time of the motor chamb er 
within the cortex! of the 2V-port model. The 
response acco r ding to the simple transfer function 
is also shown m Figure 14. It is obvious from die 
figure that die simple filling time argument for the 
hybrid combustion chamber is not valid. 


Another obser v a ti on from Figure 14 is 
that the n umerical simulation result shows 
oscillatory be ha vior in the early stages of the 
throttling event. Tins oscillatory behavior is due to 
the excitation of the acoustic modes of the cavity 
by the dis tur banc e caused by the sudden mcrease in 
the oxidizer mass flow rate. The oscillation 

acoustic vnode fundings for the same motor. 

Hus kmd of behavior is commonly 
observed in AMROC DM-OI motors 4 . Namely, a 
nonlinear d is t ur bance causes the wefl damped 
oscillation of the chamber p re s sure at the 
fongitiidina] aco u s ti c frequencies. The exerting 
disturbances in the actual tests are believed to be 
caused by: 1) a sudden blockage of the nozzle 
throat area by some disin tegr a ted fuel fractions, 2) 
the ignition of a bulk of fuel oxidizer mixture 
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captured in the pre-combustion chamber or 3) the 
low frequency oscillations of chamber pressure. 

We c om pare the predictions of our gas 
dynamic model for die acoustic frequencies of the 
chamber with the AMROC DM-01 motor test data. 
Table 7 shows the first five longitudinal 
frequencies de t er min ed by the pertuibatkn 
solutions. Since our model is isothermal, these 
frequencies are corrected by multiplying with the 

factor of . An ivc ng t j value of 1-2 is used 

for the calculations. Note that the correc t ed 
frequencies falls m the range of observed 
frequencies for every mlmtaifd acoustic mode. 

The gas d yna mi c code de v eloped to 
perform simu lati o ns on the 2V-PM model is used 
to confirm the f indin gs of the p e rtur b ati on 
solutions on the frequency response of the system. 
In the simulations, sinusoidal oxidizer mass flow 
rate inputs at various frequencies are applied as the 
forcing function and the phase and the amplitude of 
the cham ber p re ss u r e oscillations are determined. 
The agre e mcB t between the p er tur b ati on solutions 
and die n u meric al simulations for the frequency 
response of the system is also found to be 
satisfactory. 

The simulation code is also used to 
investigate the response of the gas dynamic system 
to a disturbance. The ca udiar mas flow in p ut is 

15) Tables and Figures: 


suddenly increased 10% above its operating value 
and held at its elevated level for 0.00 seconds 
before its suddenly reduced to its original value. 
The pressure time history for this simu lation is 
shown Figure 15. It is clear that die disturbance 
excites the acoustic modes of the chandler similar 
to the previous observations. These exited inodes 
decay in a moderate time scale as it was observed 
in the AMROC tests. A similar simulation is 
performed to determine the effect of nozzle 
blockage on die stability by reducing the nozzle 
throat area for a brief period of tune. The resuhs 
for the nozzle area p er flatio n case are iilfn i fji 
in nature to the resnfes of the oxidizer mass flow 
rate perturbation, hi short, a nonlinear di stu rbance 
causes the ringing of die hybrid cham ber which 
decays in time. 

Finally a simple delay (r M ) on the 
co mb u s ti on r espons e is applied to the gas dynamic 
system. 

F n =a t nG m V'*"’ (AJ21) 

Both the p ert u rba ti on solutions and also 
the numerical simulations indicate that any positive 
combustion delay have negligible effect on the 
stability nabne of the system. 





1) Solid phase kinetic times 

racier 3 sec 

Degradation mechanisms of the 
polymer 

2) Gas phase kinetic times 

< 10 3 sec 

Hydrocarbon combustion 

3) Feed system response 
times 

(Varies greatly from system to 
system) 

Response time of the feed system 

4) Evaporation times 

t^AUo^AP) 

Evaporation process of the liquid 

oxidizer 

5) Thermal lags in solid 

Tj * r/r 2 stllT 1 sec 

Thermal profile changes in the solid 
gram 

6) Boundary foyer diffusion times 

f w ac Lf U e « 10" 1 sec 
(Varies greatly form case to case) 

Turbulent boundary foyer diffusion 
processes 

3) Acoustic times (longitudinal) 

T a ac l\c « 10 3 sec 
(Varies greyly form case to case) 

Propagation of the acoustic waves 

7) Gas dynamic filling times 

Tjgg oc L/c «10 1 sec 
(Varies greatly form case to case) 

Global mass flow balance 


Table 1: The transient time scales of various phenomena in atypical hybrid rocket. Only the last four are 
modeled in this paper. 
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AMROC 

OfcnHt) 

OIF 

• 

C *o 

(ft/sec) 

* 

C c *P 

(ft/sec) 

L m 

(in) 

/ 

(Hz) 

r «2 

(msec) 

S Motor 
(10k) 

1.45 

5584 

5361 

570 

11 

48.9 

Half-Scale 

(33k) 

1.45 

5584 

5361 

1305-1430 

4.8 

117.4 

H-500 

(75k) 

1.45 

5584 

5361 

1770-2340 

4.0 

1763 

DM-01 

(250k) 

135 

5700 

5472 

2168 

2-33 

204.8 


Table 2z S ummar y of parameters used in the frequency estimations for AMROC motors. 


Test 

Nn. 

D 

(in) 

K 

(ffl/sec) 

<7. 

(Ih/ in 2 - 
sec) 

G, 

(fc/in 2 - 

sec) 

Pc 

(1*0 

OIF 

/ 

(Hz) 

*Pa,c 
( 1*0 

r «2 

(msec) 

i 

(GOX) 

435 

6J0 

0.423 

0371 

m 

385 

13 


48.8 

2 

(GOX) 

435 

■9 

0.402 

0341 

550 

389 

8 


473 

7 

(GOX) 

430 

4.7 

0339 

0.458 

m 

389 

15 

100 

44.7 

8 

(GOX) 

432 

5.1 

0365 

0.492 


387 

20 


253 

1 

(LOX) 

- 

- 

0322 

031 

435 

231 

8 

250 

663 


Table 3: Sum m ar y of p ara m e t ers used in the frequency estimations for the 1 1 inch Hybrid motor tests. 


Grain length is 102 inches. 


Test 

Na. 

D 

(in) 

L 

(in) 

(HW sec) 

G. 

(Ib/in 2 - 

sec) 

G , 

(Ih/ in 2 
ia-sec) 

Pc 

(psi) 

OIF 

/ 

(Hz) 

A**.* 

<%) 

T ui 

[msec) 

3 

■a 

102 

3.44 

038 

oszr 

745 

US 

6-10 

15 

663 

4 

43 

102 

336 

023 

0325 

325 

2.4 

10-20 

20 

47.4 

6 

■a 

102 

836 

0.78 

1.024 

750 

33 

10-15 

12 

33.7 

7 

53 

102 

3.41 

mi 

032 

335 

2.1 

6-15 

15 

733 

8 

a 

102 

3X7 

KH 

SB 

mi 

36 

2-5 

60 

191.7 

mu 

6.1 

102 

634 

038 

0392 

215 

36 

6-20 

33 

25.9 

15 

5.4 

m 

1030 

0.45 

0399 

1025 

isa 

4 

10 

83.9 

10 

63 

102 

5.93 

030 

0378 

215 

36 

10-25 

5 

363 


Table 4: Sumnaiy of p aram eters used in the frequency estimations for the NASA/MSFC 1 1 inch Hybrid 


motor tests. 
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I Hjjiriitil 

Estes* erf 

Ctrmted 

Observed 

Estimated 

AtMsticlMe 


Ficptacy 

Fmprfj 

VupBfii nlii ■ 


<B*> 

(Hz) 

(Hz) 

(l/*«) 

i* 

48.4 

53.0 

50-55 

-22 

2T 

96.7 

103.9 

100-110 

-22 

y* 

142j6 

156.1 

150-165 

-21 

4* " 

191.1 

2093 


-19 

5* 

240.7 

263 j6 

250-275 

-18 


Table 7: The csdnated smd observed longitudinal acoustic fre quencie s for AMROC DM-OI motor. 




Fipre 2. The s ch ematic of the TCG coupled system. 
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System Zero 
at (0,0) 
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20 Reaf„ 1/s ec 


Figure 3: Plot of the transfer function of the TC coupled system with no delays ( r m = r bn = 0 ). This 
is for a HTPB system with E a =1 5 kcal/mole. 



so 


Figure 4: The unstable poles of the TC coupled system in the s plane for r m = 0 , r bl2 = 38 msec. This plot 
is for a HTPB system with E a =15 kcal/mole. 
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FigiHre 5: Effect of activation energy on the oscflbtion frequency and amp lifi cation. 



Fignre Effect of blocking factor exponent on the oscillation frequency and amplification. 
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Figure 7: Effect of boundary layer delay time on the oscillation frequency for various activation energies. 
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Figure 9: The contour plot of the TCG coupled system with r m = 38 msec. Most of the transient features 
of a hybrid system can be deduced from the locations of the poles. 



FSgwre If; Chamber and feed pressure tine traces for the paraffin-based motor test, 4L-05- 
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Figure 11; Fourier transfonn of the chamber pressure for the paraffin-based motor test, 4L-05. 
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Figure 12: Comparison of the TCG coupled oscillation frequency prediction with the hybrid motor test 
data available m the literature. 
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Figure 15: The response of ihc prv-conifaeslion chamber pressure to a distixrbaioc of the oxidizer mass 
flow fate. 
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